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We discuss the impact of matter effects in the D° — D° system. We show that such effects could, 
in principle, be measured, but that they cannot be used to probe the mass difference xo or the 
lifetime difference j/d. This occurs because the mixing effects and the matter effects decouple at 
short times. We also comment briefly on the B systems. 
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I. INTRODUCTION 



o 
o 
o 

£S| ■ The D° — D° system is of great current interest. Until recently, the best probes on the mixing resulted only in 
upper bounds on x 2 D +yp of order 1CU 2 This changed with the recent results from CLEO B and FOCUS ||. The 
current situation is beautifully summarized in reference Q|: there is evidence for yjj ~ 1CP and for a large strong 
phase, corresponding to large SU(3) breaking effects. The problem with the strong, final-state interaction phases is 
ON , that they are not known. They are fit in the same experiments that look for yo- It would be interesting if matter 
effects could be used to shed some light on this issue. The idea arises from the kaon system, where the phases due to 
matter effects are measured experimentally. Thus, we would be dealing with a known CP-even phase. 

There is one important difference between the kaon and D systems. For a beam of neutral kaons in vacuum, the 
t^J- \ depletion from the beam is controlled by two exponentials, exp (— Ts t) and exp (— Tl i), corresponding to the short- 
lived (Ks) and long-lived (Kl) components. Here |AT| ~ 2T and the two exponential fallouts are clearly separated. 
As a result, we may wait for the Ks component to decay away and use the matter effects on Kl to regenerate K$; 
a phenomenon known as 'regeneration'. In the D (and the Bd) system the situation is very different: |Ar| <C T and 
the leading behavior in vacuum is given by exp(— Tt), with the hyperbolic sine and cosine of ATt acting as small 
perturbations. In this system, the depletion from the beam is controlled by T. Because we cannot resolve the two 
exponentials, the classical 'regeneration' experiments cannot be carried out and a new analysis is required. This is 
i-£h ' the question we address here. 

In section [n| we review the general features of propagation of neutral mesons systems in matter. In section III 



we 

O ' use an extrapolation from known kaon results to argue that these effects might, in principle, be sought. Unfortunately, 
we will also show that, even if visible, the matter effects cannot be used to probe xd or yo- This is shown to be a 
' generic feature of the 'small-time' approximation in the time-evolution. We draw our conclusions in section [V, where 
\ we remark briefly on the and B s systems. 

•rH , 

X: 

II. PROPAGATION IN MATTER 

Let us consider a generic neutral meson system P a — P°. Assuming CPT-conservation in vacuum, this system is 
characterized by two complex eigenvalues, fi a = m a — and /j,b = — iTb/2, and by a mixing parameter, q/p. 

It is convenient to define the average and the difference of the eigenvalues, 

Ma + M6 i-n 

ii = = m 1 . 

p 2 2 ' 

i 

A/i = fi a - [ib = Am - 2^ r - (!) 

These quantities describe the system in vacuum and its time evolution is well known [pi . 

The new effects resulting from the interaction with matter can be described by the elastic forward scattering 
amplitudes of P° and P°, which we denote by / and /, respectively. These enter the effective Hamiltonian through 
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where N is the density of scattering centers in the medium. As before, we define the average and the difference of 
these parameters as 



Xa 



x + x 



Ax = x - X 



2ttN 



-A/, 



(3) 



will: Af / /. 

In the presence of matter, the time evolution of meson states with a given initial flavor may be written as 



e *(M + Xav)t |p0^ = 



Ay 

cos (fi t/2) - %-± sin (fi t/2) 
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Ay 

cos (fi t/2) + i-^ sin (fi t/2) 



(4) 



where fi = ^(A^) 2 + (A x ) 2 @. The evolution in vacuum is reproduced by setting Xav and Ax to zero. 

As mentioned in the introduction, in the D (and in the Bj) system the depletion from the beam in vacuum is 
controlled by T. This one exponential is then modulated by trigonometric functions of A/jt = -j£ r, where t = Tt is 
the time measured in lifetime units. Therefore, one introduces 



x — ly. 



(5) 



In most cases of interest the matter effects satisfy T > |lmxav| and A/i > |AyL and the depletion is still chiefly 
determined by T. To describe the trigonometric functions in matter we need Eq. (|5|) and a new complex parameter. 
This may be choosen to be A\/T, which determines the depletion (statistics available) by the time the effects of Ax 
become important. However, since we know the current experimental reach on Afi, it is easiest to use it to access the 
experimental reach on Ax, through^ 



2r 



Ax 
Afi' 



(6) 



Our figure of merit will be \Ax/A[i\. 

Henceforth we will refer to those effects due to A/i as 'mixing effects', and to those effects due to Ax as 'matter 
effects'. The fundamental question is: when do the matter effects and the mixing effects couple to each other? The 
answer is easily found by looking at Eqs. (^). They can couple in two ways. The first and simplest way occurs when 
the production, propagation and decay occur all in the same medium and the decay is flavor tagging. In that case, 
the matter and mixing effects couple mainly through fi. The second way occurs when both flavors can decay into 
the same final state, or when the system crosses several media. For example, for events occuring completely in the 
medium, with the decays going into CP-eigenstates. Or, when the production and the initial part of the evolution 
occur in matter, with the final part of the evolution and its decay occuring in vacuum. In these cases, besides fi, 
there can also be effects proportional to Afj, Ax, which arise from the interference between the coefficients of the sine 
terms multiplying and of the sine terms multiplying \P°). 

III. MATTER EFFECTS IN THE D° -D° SYSTEM 

A. Interplay between mixing and matter effects 

We will now show that, even when visible, the matter effects cannot be used to probe Xd ot Ud- The main reason 
is that the complex trigonometric functions are controled by (x — iy)r and, with current technology, one can only 
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follow r up to 5 or 10, corresponding to 5 to 10 average lifetimes in the D proper frame. This is mostly due to the 
statistics of perfectly reconstructed events. It is obvious that, given unlimmited event samples, one could follow the 
full decay curve including times of order 1/ud- We mention this because the 'no-go theorem' we are about to present 
hinges on the experimental constraint that one can only follow the decay curves for time t <c 1/ud- Under these 
circumstances, we can expand the trigonometric functions on the right-hand side of Eqs. (^) as 

c„ s ( n m ± if ».» (a ( /2) - 1 ± if t - e T ,a x <AE£±mi ,» + 0{l% (7) 

and 

- ^ sin (fl t/2) *~i^t + iA, (^) a + W ,a + Qn (g) 

for the right-sign and wrong-sign transitions, respectively. 

The old results on semileptonic wrong-sign decays jll] measure the magnitude squared of Eq. (|^) , given to leading 
order by 

4r 2 4 ( ' 

Let us now consider a final state / to which both P° and P° can decay. We denote the decay amplitudes by Af and 
A/, respectively. Such a decay has a term linear in t arising from the interference between Eqs. (0) and 



Im 



A^A f A) 



(10) 



CLEO (2| has looked at / = K + n which is a Cabibbo-allowed decay for D° and a Doubly-Cabibbo-suppressed decay 
for D°. FOCUS § has looked at the CP-eigenstate / = K + K~. 

Could one look at matter effects with semileptonic decays? In principle yes. The time evolution of the right-sign 
semileptonic decays is proportional to the magnitude squared of Eq. (^), which has a term linear in t, 

Tlm(y) r. (11) 

Even if |Ax| is an order of magnitude smaller than \Afi\, this term should still be easier to detect than the quadratic 
term in Eq. (||). This idea will be explained in detail below. 

For the moment we wish to return to our fundamental question: can matter effects be used to give a new handle 
on A/i? The answer is no! First consider right-sign semileptonic decays and magnitude square the right-hand side of 
Eq. (pj). We see that only at order r 3 do we find effects proportional to A/i A%. Even at order t 2 , the effects decouple 
as Re(A/i) 2 and Re(Ax) 2 . Next consider wrong-sign semileptonic decays and magnitude square the right-hand side 
of Eq. (||). This case is even worse, since the effects proportional to A/j, A\ only arise at order r 4 . Now consider again 
a final state / to which both D° and D° can decay. In the presence of matter Eq. ( |Io|) gets changed into 



Im 



.iA lA) {t + if-e 



(12) 



In this case there is an effect proportional to A/z A%* at order t 2 . But, if \A\\ is an order of magnitude smaller than 
|A/x|, then this effect will be an order of magnitude smaller than that in Eq. Notice that, in particular, there is 
no effect proportional to fit. 

In this discussion we have implicitly considered only events which occur totally in matter. However, the matter- 
vacuum transition can only promote the same type of interference effects we have already discussed p2[ . 

B. Observing matter effects 

We have already argued that matter effects cannot be used as a handle on and yo. We have yet to prove that 
they are observable at all, even in principle. Let us consider the decay D° — > K~l + i>i. Using Eqs. (Q) and (0) we find 

T[D°{t) -^K-l+n] «e- rt e 2lm ^*|A | 2 (l + ImA X t), (13) 
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where we have denoted A a = A{D° -» if-f+iyj)^ We now need to estimate the linear term. 

However, since the elastic forward scattering amplitudes for D and D° (/ and /) are not known, we have to estimate 
them somehow. We rec all that the matter effects in the kaon system arise because K° interacts quasi-elastically with 
the nucleons, while K° suffers inelastic interactions, such as 

W + (p+ n) ^tt°+ (£+ £°) 
sd + (uud,udd) — » dd + (uus,uds). (14) 

The counterpart in the D system is 

D° + (p+,n) ^vr + (E+ £°) 
cw + (uud,udd) — > + (udc,ddc). (15) 

The kinematics are different because m(D°) w 1865MeV > m(ivT ) « 498MeV and m(S+) w 2454MeV > m(E+) « 
1189MeV. However, on the one hand we will be interested in the high momentum range (larger than 25GeV), where 
the difference becomes less important. And, on the other hand, we are not interested in exact results but, rather, 
wish to learn whether the matter effects are in principle observable. Therefore, we will simply scale the results from 
the kaon system to the D system. 

We will use an empirical scaling law determined by Gsponer and collaborators from their measurements in C, Al, 
Cu, Sn, and Pb for kaon momenta between 30 to 150 GeV. They found that 

/ A x 0.758 386 

|A/| = 1.13fm — 7-B-t) > (16) 



L4J 

in 




gmor 1 / VGeVc- 1 

where A the atomic number of the material. This result exhibits a power law momentum dependence in accordance 
with Regge theory [jl5|, which also predicts that arg A/ should be constant and given by —(1 + 0.386)7r/2 = — 0.6937T 
The power-law approximation is fairly good down to a few GeV/c momentum, where low-energy resonances set 
(I. We neglect these effects since we will be concentrating on the high momentum range. We will also need the 
imaginary part of / + /. Again we rely on Gsponer and collaborators, who found |18 1^] 

/ A \ °- 840 

I m (/ + /) = l. 8 95 fm (_ pT ) (g^) ,17) 

in the kaon sector. 

We start by noticing that the density of scattering centers is a medium is given by 

N=^f, (18) 

where Na is Avogadro's number, p is the density of the material, and A is its atomic number. Using Eqs. (^) and 
(|l6|), we obtain 



A XD 1 / 2ttN aP 



. Af 

/ a \ -0.242 / x -0.242 

= -5.6xl0-vW^_) (-JL-) . (19) 

VgmorV \gcm-3) VGeVc- 1 / v ' 

This result should be compared with that obtained for the neutral kaons, where —5.6 x 1 0~ 5 is substituted by 0.09. 
The sign difference arises from the fact that the inelastic interactions occur for D°, not D°; for example, Eq. ( |T5| ) 
involves D°, while Eq. ([L4|) involves K°. Since we are assuming that the forward scattering amplitudes coincide, the 
difference in magnitudes arises from 



2 To illustrate this effect, we have simplified the discussion by considering a 'thought-experiment' with the events taking place 
(production, evolution, and decay) inside the material. For these purposes, a muon in the final state would be most readily 
observed. 
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|Ax£)|/r D m K T K 498MeV 0.415 x 10~ 12 s ^ t , 



|AxK|/r K m^rz, 1865McV2 x 0.893 x 10" 



= 6.2 x 10~ 4 . (20) 



The results for carbon and tungsten are listed in table I for a variety of D momenta chosen within a range accessible 
to FOCUS JlTy - Tungsten is a preferable material since it leads to values which are roughly 4.4 times larger than 
those obtained in carbon. For tungsten and momenta above lOOGeVc -1 the effects are of order 0.2%. However, the 
current experiments already probe A/io/^d of order 1%. Therefore, in principle, the linear term in Eq. ( |l3| ) would 
be detected with an experiment of this type. 

Of course, we still need to show that the absorption term exp(2Imx av i) does not affect greatly the number of 
events. Using Eq. (|l7|), we find 

2 lmX - D = -9.4 x 10- 5 ( A T ) ~°' 16 ° ( -J>) ( PD -) . (21) 
T D VgmorV \gcm-V VGeVc" 1 / V ' 

For tungsten and a momentum of lOOGeVc" 1 , the ratio becomes 0.08. Looking back at the right-hand side of Eq. (fl3|), 
this result means that the absorption exponential, exp (21m % av t) , does not compete with the decay exponential, 
exp(-rt). 

Although we have been detailed about our estimates, our objective is not to propose a specific experiment, but 
rather to emphasize that one is not, in principle, out of reach. Therefore, we will not address the obvious list of 
questions: i) How does one identify the primary and secondary vertices?; ii) How does one instrument the material?; 
iii) What is the ideal momentum range?; etc. We stress that the results presented here are based on the expectation 
that we can extrapolate from the kaon to the D system. They should be taken as estimates and not as reliable 
predictions. In particular, the presence of resonances could enhance or suppress these estimates. 

IV. CONCLUSIONS 

The situation in the Bd system is similar to the one in the D system in that Ar <C T. Using the same naive 



high- momenta, scaling law used in Eq. (20), we obtain 



\Axb\/T b _ m K T K 
\Axk\/Tk m B T B 



= 8.2 x 10~ 4 . (22) 



However, in the Bd system |A/i| ~ Xd ~ 0.7, meaning that |A%| cannot compete with V nor with |A/i|. 

The situation in the B s is very different. From the point of view of the lifetime difference it is actually closer in 
spirit to the kaon system pb]| . Indeed, in the B s system |Ar|/T ~ 0.15 |2l}| , and one may talk about two exponentials. 
However, B s contains no quarks of the first family and, hence, it would seem inappropriate to estimate the matter 
effects as we have done above. 

In conclusion, we have analyzed the possibility of detecting matter effects in the D° — D° system. We conclude 
that these effects are, in principle, accessible. Unfortunately, when the time dependence is followed only for a few 
lifetimes, the matter effects and the mixing effects decouple. As a result, matter effects cannot be used to provide a 
new handle on xd or yr>. 
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in tungsten 


GeVc-i 




r D 


25 




0.024 




0.106 


50 




0.032 




0.138 


100 




0.042 




0.182 


150 




0.048 




0.212 


200 




0.054 




0.236 


250 




0.058 




0.258 



TABLE I. Values of Axd/^d for carbon and tungsten, as a function of the D momentum. 
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